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ABSTRACT 

We consider the degree to which "21 cm tomography" of the high-redshift Universe can distinguish different 
ionization histories. Using a new analytic model for the size distribution of H II regions that associates these 
ionized bubbles with large-scale galaxy overdensities, we compute the angular power spectrum and other statis- 
tical properties of the 21 cm brightness temperature during reionization. We show that the H II regions imprint 
features on the power spectrum that allow us to separate histories with discrete bubbles from those with partial 
uniform ionization (by, for example, X-rays). We also show that "double reionization" scenarios will modify 
the morphology of the bubbles in ways that depend on the mechanism through which the first generation of 
sources shuts off. If, for example, the transition occurs globally at a fixed redshift, the first generation imprints 
a persistent feature on the 21 cm power spectrum. Finally, we compare our model to one in which voids are 
ionized first. While the power spectra of these two models are qualitatively similar, we show that the underly- 
ing distributions of neutral hydrogen differ dramatically and suggest that other statistical tests can distinguish 
them. The next generation of low-frequency radio telescopes will have the sensitivity to distinguish all of these 
models and strongly constrain the history and morphology of reionization. 
Subject headings: cosmology: theory - intergalactic medium - diffuse radiation 



1. INTRODUCTION 

Recently, a great deal of effort - both observational and 
theoretical - has been focused on understanding the reioniza- 
tion of hydrogen in the intergalactic medium (IGM) at z > 6. 
Several independent observational techniques offer comple- 
mentary views of this landmark event. The most straightfor- 
ward method is to seek extended troughs of complete Lya 
absorpt ion in the spectra of high -redshift quasars. Evidence 
for this Gunn & Peterson ( 1965) effect has been found near 
6.5 (Be cker et alJ2001llWhite et al.l2003l) . Unfortunately, 
the optical depth of a fully neutral medium is so high that 
current measurements only require a mean neutra l fraction 
xh > 10~ 3 , and even this limit has been disputed (Songaila 
2004). Studies of the rapidly growing ionizing background 
(Fan et al. 2002 ) and the Stromgren spheres surrounding the 
quasars ( Wvith e & Loebl 2004) suggest larger neutral frac- 
tions but depend on detailed modeling. A second constraint 
arises because free electrons scatter cosmic microwave back- 
ground (CMB) photons, washing out intri nsic anisot ropies 
and generating a polarization signal (Zaldarriaga 1997). This 
provides an integral constraint on the reionization history; re- 
cent measurements of the C MB imply that xh mus t have been 
small at z > 14 (iKogut et al.H200l LSpergel et alJl200l . Fi- 
nally, the relatively high temperature of the Lya forest at 
z ~ 2-4 suggests an order unity change in xh at z < 10 
( iTheuns et al.ll2002l iHui & Ha iman 2003), although this ar- 
gument depen ds on the c haracteristics of He II reionization 
(e.g.. LSokasian et all2002l) . 

Taken together, these three sets of observations suggest a 
complex ionization history extending over a large redshift in- 
terval (Az ~ 10). This is inconsistent with a gen eric picture 
of fast reionization (e.g.. iBarkana & Loebl l2001. and refer- 
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ences therein). The results seem to indicate strong evolu- 
tion in the sources responsible for reionization, and a precise 
measurement of the ionization h istory would strongly con- 
strain early structure formatio n (IWvithe & Lqebl|2" 003l ICenl 
I20031 lHaiman & Holderll2003t ISokasian et alJl2003bl Un- 
fortunately, extracting such detailed information requires new 
observational techniques. One of the most exciting pos- 
sibiliti es is "21 cm tomograp hy" of the high-red shift IGM 
dScott&Reeslll990t iKumar etal.ll 19951 iMadau et alJl!997l) . 
in which one maps the distribution of neutral hydrogen on 
large scales through its hyperfine transition. By probing a spe- 
cific spectral line, 21 cm surveys measure the time history of 
reionization, and unlike the Lya forest they do not suffer from 
saturation problems. In fact, the signal is sufficiently weak 
that foregrounds from the Galaxy and a wide variety of ex- 
tragalactic obje cts pose su bstantial challenges to these experi- 
ments (Oh & Mackl2003tlDi Matteo et alJ2004l) . Fortunately, 
the known foregrounds all have smooth continuum spectra, 
which should allow frequency cleaning to high accuracy be- 
cause the 21 cm signal has spectral struct ure on small scales 
dZaldarriaga et al.ll2004t hereafter ZFH04; iMorales & Hewitt! 
l2003UCoorav & Furlanettol2004l) . 

However, predicting the 21 cm signals from reionization 
has proven difficult. One method is to use numerical si mula- 
tions of reionization dCiardi & Madaull2003l iFurlanetto et alJ 
2004a; Gnedin & Shaver 2003), but computational costs have 
so far limited the simulations to subtend (at best) a few reso- 
lution elements of the 21 cm maps. Analytic models of reion- 
ization can extend to larger scales but require some way to 
describe the complexities of structure formation and radia- 
tive transfer. The simplest approach, in which we assign each 
galaxy its own H II region, does a poor j ob of m atching the 
large ionized bu bbles found in simulations ( Ciard i et alJ2003l 
Sokas ian et alJl2003al) . We have recently developed a model 
that associates H II regions with large-scale fluctuations in the 
density field and reproduces the major features seen in sim- 
ulations (Furlanetto et al. 2004b, hereafter Paper I). This al- 
lows us, for the first time, to make quantitative predictions 
about the 21 cm signal at reionization. With high signal-to- 
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noise 21 cm maps, we can directly measure the distribution of 
H II regions in order to test our model. However, the signals 
are sufficiently weak that such maps will be difficult to make; 
fortunately, statistical measurements of the data also contain a 
great deal of information about reionization (ZFH04). In this 
paper, we use the model of Paper I to predict the 2 1 cm an- 
gular power spectrum for several reionization scenarios, with 
an emphasis on how the measurements help distinguish the 
crucial features of these different pictures. In {j5]we briefly 
review our model for reionization. We then show how 21 cm 
measurements can distinguish different reionization histories 
in ^5] and how they can distinguish different models of reion- 
ization in |4] We conclude in ^5] 

In our numerical calculations, we assume a ACDM 
cosmology with tt m = 0.3, fl A = 0.7, Sl b = 0.046, H = 
100/zkms" 1 Mpc" 1 (with h = 0.7), n = 1, and g« = 0.9, consis - 
tent with the most recent measurements ( Sn ergel et alJ2003l) . 

2. A MODEL FOR REIONIZATION 

Recent numerical simulations (e.g.. ISokasian et a"" l2003a) 
show that reionization proceeds "inside-out" from high den- 
sity clusters of sources to voids, at least when the sources 
resemble star-forming galaxies (e.g., Springel & Hernquist 
2003). We therefore associate H II regions with large-scale 
overdensities. We assume that a galaxy of mass m ga i can ion- 
ize a mass (mgai. where £ is a constant that depends on the 
efficiency of ionizing photon production, the escape fraction, 
the star formation efficiency, and the number of recombina- 
tions. Values of £ ~ 10—40 are reasonable for normal star 
formation, but very massiv e stars can increase th e efficiency 
by an order of magnitude (Bromm et al. 2001b). The crite- 
rion for a region to be ionized by the galaxies contained in- 
side it is then / co n > £ _1 , where / co ii is the fraction of mass 
bound to halos above some m m \„. We will normally assume 
that this minimum mass corresponds to a virial temperature 
of 10 4 K, at which point hydrogen line coolin g becomes effi- 
cient. In the extended Press-Schechter model ( L acev & Cold 
1993), this places a condition on the mean overdensity within 
a region of mass m, 
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where K(Q = err'(l -C" 1 ), cr 2 (m) is the variance of density 
fluctuations on the scale m, <r^ in = <r 2 (m m m) and 5 c (z) is the 
critical density for collapse. 

Paper I showed how to construct the mass function of H II 
regions from 5 r in a n anal ogous way to the halo mass function 
JPress & Schechterl ll 9*741 Bond et al.1 f 19911) . The barrier in 
equation Q is well approximated by a linear function, 5 X « 
B(m,z) = Bo+Bia 2 ( m). In that case the mass function has an 
analytic expression (Sheth 1998): 
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where p is the mean density of the universe. Equation (0 
gives the comoving number density of H II regions with 
masses in the range (m,m + dm). The crucial difference be- 
tween this formula and the standard Press-Schechter mass 
function arises from the fact that 5 X is a (decreasing) func- 
tion of m. The barrier is more difficult to cross as one goes to 
smaller scales, which gives the bubbles a characteristic size. 
In contrast, the barrier used in constructing the halo mass 
function is independent of mass, yielding the usual power law 
behavior at small masses. The characteristic scale depends 
primarily oni/f. 



Paper I also showed how to construct the power spectrum of 
the quantity ip = xh(1 +5). We explicitly computed the correla- 
tion function of the H II regions £ xx in terms of the mass func- 
tion in equation (0, and we computed the density correlation 
function £gg and cross correlation between de nsity and ioniza- 
tion f raction £,5 using the halo model (e.g.. lCoorav & Sheth 
120021) . We combined these correlation functions to get the 
correlation function of tp, 
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ZFH04 showed how to convert this to the three-dimensional 
power spectrum P^ and then to the (observable) angular 
power spectrum of the 21 cm brightness temperature 5Tb, es- 
sentially by doing a Fourier transform. Here 

JT 00 , {T s -Tcmb\ ( £lhh 2s 
5Tb ~ 23 ip x 
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with Ts the spin temperature of the gas. In all of the calcula- 
tions here, we assume Ts ^> 7cmb, which should be reasonable 
even relatively early in reionization (see §2 of ZFH04). We 
further assume perfect frequency resolution in computing the 
angular power spectrum; our results therefore do not depend 
on the experimental setup. The effects of finite bandwidth 
are described in ZFH04 and Paper I. We refer the reader to 
these papers for more details on our approach. Note that these 
methods to construct the angular power spectrum of 5Tb can 
be applied to any model of the H II regions. 

3. COMPLEX REIONIZATION 

Paper I examined the power spectra of standard reionization 
scenarios in which all the ionizations are confined to discrete 
H II regions. We showed that the power spectrum evolves 
rapidly throughout reionization, allowing one to reconstruct 
the ionization history in such a simple scenario. Here we ex- 
tend our model in order to consider the effects of more com- 
plicated reionization processes. 

3. 1 . A Partially Ionized IGM 

We first consider how well 21 cm observations can dis- 
tinguish a model with all ionizations inside of H II regions 
from one with a uniform com ponent, such as may h appen 
through ionization by X-rays JOhl 1200 It [Venkatesan et alJ 
120011: iRicotti & OstrikeJl2003HMadau et alJl2003l) or decav- 
ing particles ( Sciamal [l982t iHansen & Haimanl 120041) . Sup- 
pose that the IGM has a uniform ionization fraction x u (z). 
On top of this uniform level are some variations in the neu- 
tral fraction owing to isolated H II regions. In this case, the 
condition / co u > £ (1 — x u ) replaces the barrier in our model: 
each galaxy can produce a larger ionized bubble with the same 
number of ionizing photons. However, rather than varying 
from zero to unity, we have < xh < (1 — x u ). This damps 
the fluctuations from the bubbles, requiring the replacements 
-> (l--x„) 2 £« and £ xS -> (1 -x u )£, x s- Otherwise our for- 
malism is unchanged. 

Figure [l] contrasts such a case with a "standard" history. 
The thin curves assume £ = 40 and x u = 0, while the thick 
curves have £ = 12 and set x u (z) to force the total neutral frac- 
tion to match in the two cases. The uniform component is thus 
responsible for 70% of the ionizations. When Xh « 1, the two 
curves are nearly identical because the bubble feature is in- 
trinsically weak. However, while the bubbles rapidly imprint 
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FIG. 1. — The 21 cm fluctuation amplitude as a function of redshift. The 
thick curves have £ = 12 plus a uniformly ionized component, while the thin 
curves have £ = 40. The ionization fraction xh (z) is identical in the two mod- 
els. The curves have: xh = 0.96 (dotted), xh = 0.8 (dashed) and J# = 0.26 
(solid). 



a feature on the power spectrum if x u = 0, uniform ionization 
strongly damps the power from the H II regions. As a result, 
in such a scenario the 21 cm fluctuations trace those of the 
density (normalized by 1 - x u ), at least until reionization is 
almost complete. Consequently, the amplitude of 5T/, is sup- 
pressed by 25-50% on scales near the characteristic bubble 
size, even when xh is large. We find that, unless reionization 
is near completion, the bubbles need to provide a majority of 
the ionizations in order for them to have an appreciable effect 
on the power spectrum. This provides a clear test for whether 
Stromgren spheres from ultraviolet photons or smoother ion- 
ization from some other source dominate the morphology of 
reionization, even if the survey is unable to make high signal- 
to-noise maps. 

3.2. Double Reionization 

Recently, a number of theoretical models for reionization 
have attempted to reconcile the CMB and quasar data by pos- 
tulating an early generation of sources with high ionizing ef- 
ficiency (most often because they contain massive, metal free 
stars) along with a self-regulation mechanism that switches 
to normal star format i on with a l ower ionizing efficiency 
jWvithe&Loebll2003t ICmll2003t lHaiman & HoldeJ 120031 
ISokasian et al.ll2003bl) . Such scenarios can lead to a plateau 
in the ionized fraction or even "double" reionization, in which 
ionized phases bracket a mostly neutral period. In this section, 
we explore the signatures such scenarios imprint on the 21 cm 
power spectrum through some simple toy models. 

We consider an early generation of sources with (Ci ,fWmin,i) 
and a later one with (C2,fWmin,2)- In the initial phase only the 
first generation is important, so the ionization simply follows 
the model of Paper I. We then end this generation according 
to one of two conditions. The first choice shuts down these 
sources once the universe reaches some uniform level of ion- 
ization at redshift z\. For example, the early phase coul d be 
caused by X-rays or decaying particles, as described in 33.11 
A better- motivated scenario could be the following. Sup- 
pose that the first generation consisted of massive, metal-free 



stars; the natural self-regulation condition halts the forma- 
tion of these stars when the metallicity in co l lapsed objects 
passe s some thresh old (|Bromm et alJ 12001 at iMackev et alJ 
2003; Yoshidaetal. 2003). These sources will stop forming 
soon after they first appear in any given region (or in other 
words soon after the region is ionized for the first time) but 
they will continue to form in the remaining neutral regions. 
Thus the first generation will persist until the entire universe 
has been ionized. Our toy model makes the simplification that 
the first generation was able to keep the universe ionized un- 
til the transition occurs; we implicitly neglect recombinations 
within regions that were ionized before z\ ■ 

While the second generation is still rare, recombination in 
the IGM will be the dominant effect. Such an era is obvi- 
ously easy to identify with 21 cm tomography because it is 
accompanied by a corresponding increase in xh and hence in 
the signal tFurlanetto et a i1l2004ah . If the recombination rate 
were uniform throughout the IGM, we would have the power 
spectrum P^, = xjjPgg. Because the recombination rate is a 
function of the density, we also get fluctuations in the neu- 
tral fraction that amplify ST/,. In the limit that the elapsed 
time is smaller than the recombination time, we would have 
Pxx ~ [1 -x~i(z)/x~i(zi)] 2 Ps5, where Xi is the mean ionized frac- 
tion. Unfortunately, Xh(S,z) from recombinations is a nonlin- 
ear function, so the detailed evolution must await numerical 
simulations. Given the weak large-scale fluctuations present 
at high-redshifts, we do not expect recombinations to amplify 
the power by more than a few tens of percent, and it should not 
imprint prominent features into the power spectrum. (Note 
that allowing extra recombinations in regions that were ion- 
ized before z\ would increase this amplification factor.) In 
our for malism , we can implement this scenario using the re- 
sults of 33.11 with the extra wrinkle that x u can decrease with 
cosmic time. 

Another prescription shuts off the first generation before 
those sources can fully ionize the universe. At zi, the uni- 
verse has developed a patchwork of H II regions but they have 
not yet overlapped completely. This choice is appropriate if 
the self-regulation mechanism is not local. For example, the 
buildup of an ultraviolet background in the Lyman- Werner 
bands could halt H? cooling and shut off star formation in 
small halos ( Haiman et aQ 119971) . or an X-ray background 
could heat the universe and raise the Jeans mass. In this 
model, the first generation imprints a set of ionized bubbles, 
within which most of the second generation sources grow (be- 
cause both appear in the same overdense regions). The bub- 
bles grow only slowly until the total number of ionizations 
from the second generation become comparable to that of the 
first; after this point the evolution approaches the normal be- 
havior. If we neglect recombinations within H II regions, the 
barrier can be set up as follows. Consider a region of mass 
m at Z2 < Zi! as in Paper I we need the condition for sources 
inside this region to ionize it. The number of ionizations is 
simply the sum of those from both generations. Thus we find 
that the excursion set barrier S x (m) will be the solution of 



S c (zi)-S x (m) 
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(5) 



Here the complementary error functions are the fraction of 
collapsed gas above the mass thresholds at the two redshifts. 
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FIG. 2. — The 21 cm power spectrum in three different models of "double 
reionization." The solid curves assume that the universe has x„ = 0.5, the 
long-dashed curves assume zi = 18 with fi = 500, and the short-dashed curves 
assume z\ = 18 with c^i = 105 and m mm ten times smaller than our default 
value. In all cases, £2 = 12. The thin curves are for z = 16 {xa = 0.54) and the 
thick curves are for z = 12 (xh = 0.73). 



This is a more complicated expression than in Paper I, but 
8 X is still approximately linear in a 2 . With the new barrier, 
the formalism from Paper I carries over without further mod- 
ification. We could crudely incorporate recombinations into 
this model by decreasing (j as redshift decreases. This is not 
exactly correct, because it would cause the bubbles to shrink 
rather than recombining at all radii. However, this may not be 
a serious problem because, according to our prescription, both 
types of sources will primarily form within the same highly- 
biased regions. Even a relatively small number of second 
generation ionizing sources could be enough to halt recom- 
binations within the bubble unless the recombination time is 
much smaller than the Hubble time. In any case, the main 
feature of this model is that the first generation bubbles im- 
print a particular scale on the power spectrum. If recombina- 
tion is significant (but not complete) and the second genera- 
tion sources reside in partially-ionized bubbles, their own H II 
regions will grow faster than normal and will quickly reach 
this scale, at which point their expansion will slow. If, on the 
other hand, recombination is nearly complete in these regions, 
all evidence of the first generation disappears and the single 
generation description of Paper I would apply. 

Figure |2 contrasts the power spectra of three scenarios for 
double reionization. In the solid curves, the first sources uni- 
formly ionize the IGM to x u = 0.5 (or, if we neglect ampli- 
fication due to recombinations, the IGM has recombined to 
this point after a phase of full ionization). In the other mod- 
els, we end the first generation at z\ = 18 when xh = 0.5; the 
dotted and dashed curves correspond to two different sets of 
parameters for the first generation of sources. In all cases the 
second generation has £2 = 12 and ionizes ~ 4% (z = 16) and 
~ 23% (z = 12) of the IGM. Note that we have ignored re- 
combinations after z\ in order to isolate the effects of different 
prescriptions. The most obvious feature of the second type of 
transition is that the bubble size remains approximately con- 
stant with redshift. This contrasts sharply with scenarios in 
which a single set of sources is responsible for reionization, 



for which the bubble scale changes rapidly over Az ~ 1 (Paper 
I). We suggest that this plateau in the bubble scale is a clear in- 
dication that the transition between generations occurs rapidly 
and in a non-local manner. It is, however, difficult to distin- 
guish between the two scenarios containing relic bubbles. As 
described in Paper I, this is because the scale of the H II re- 
gions depends only on the behavior of the large scale density 
fluctuations and is fixed almost entirely by xh- In the uniform 
model, the amplitude of 5Tb is significantly smal ler and the 
bubble features are suppressed, just like in 33.11 Note also 
that in this case the bubble feature grows rapidly, because the 
H II regions expand into a pre-ionized medium. This stands in 
stark contrast to the dashed curves. While including recom- 
bination could increase the overall amplitude of the uniform 
model by a small amount, it would not affect the bubble fea- 
ture. 

4. OUTSIDE-IN REIONIZATION 

In Paper I and in |3] we have examined how the parameters 
of the ionizing sources change the 21 cm power spectrum; 
however, we have always worked within the model described 
in |2] Many other models for reionization exist, and we now 
consider whether the 21 cm signal can discriminate between 
them. The most popular mod el assumes that the recombi - 
nation rate controls ionization ( Miralda-Escude et al. 2000). 
Thus gas elements below a fixed density threshold are ion- 
ized first, with the threshold increasing as reionization pro- 
gresses. We refer to this as "outside-in" reionization because 
voids are ionized before the dense regions, where galaxies 
sit. A simple way to describe such a process is to assume 
that regions below a fixed mean overdensity S v < are ion- 
ized. We can then construct the mass function of the ion- 
ized bubbles throug h the excursion set form alism; the result 
is identical to the Press & Schechtejl i 19741) halo mass func- 
tion except that \S V \ replaces the collapse threshold. 4 We fix 
6 V by requiring the fraction of mass in ionized regions with 
m > (m m { D (i.e., larger than an H II region around the smallest 
allowed galaxy) to equal C/coii- This prescription is un doubt- 
edly an oversimplification of the Miralda ~Escude et alJ ( 2000) 
picture, because it completely neglects the source locations, 
the cosmic web, etc. But it nonetheless serves to illustrate the 
importance of our basic assumptions. The only additional in- 
gredient we need to construct the power spectrum is the bias 
of ionized regions. With the above prescription, we have 

b v (m,z) = -\l+ v/ \. (6) 

This is the same expression as for the halo bias in the Press- 
Schechter model ( Mo & Whita [l99a) . except that the bias is 
negative because the number density of voids increases in un- 
derdense regions. 5 

We compare the outside-in angular power spectrum to the 
inside-out model of Paper I in Figure|3] We have normalized 
the two models to have identical x~h(z)- We see that, until rel- 
atively late in reionization, the fluctuation amplitude and the 
characteristic bubble scale are similar in the two models. This 

4 In principle, we should incl ude the "void-in-cloud" process described 
by Sheth & van Wevgaert 1 2003 ) in order to remove those voids that have 
been consumed by collapsed objects. However, the regime in which we are 
interested has |<5„ <C S c , so this process turns out to have no effect on our 
results. 

5 With this prescription, we can have the cross-correlation between a halo 
and a void < -1, which is unphysical. In these situations we set = - I; 
fortunately, this has a negligible effect on our results. 
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FIG. 3. — The 21cm power spectrum in the inside-out model of Paper I 
(thick lines) and the outside-in model in which voids are ionized first (thin 
lines). In both cases, f = 40. The curves are: xh = 0.96 (z = 18, dotted), 
x H = 0.8 (z = 15, dashed), and x H = 0.26 (z = 12, solid). 



is because both models associate H II regions with large-scale 
fluctuations in the density field. The only difference is that 
the void model barrier S v is independent of scale but 6 Z in- 
creases nearly linearly with a 1 in the inside-out model. As 
a result, the void model has somewhat more ionized gas in 
smaller regions and a somewhat smaller characteristic scale; 
however, these differences may not be robust given the ap- 
proximate methods used to construct the size distribution. On 
the other hand, the similarity shows explicitly that, regard- 
less of the detailed model, the characteristic bubble size is 
large (and in the ideal range for 21 cm tomography) so long 
as the sizes of H II regions are determined by fluctuations in 
the large-scale density field. A more fundamental difference 
between the two models is that the void model has signifi- 
cantly more large-scale power, at least until shot noise from 
the bubbles begins to dominate at z ~ 13. 6 By associating 
ionized regions with underdensities, the void model amplifies 
large-scale power. In contrast, inside-out reionization washes 
out large scale power by ionizing overdensities first. 

Moreover, as emphasized in Paper I, the power spectrum 
does not uniquely determine the statistics of the 5Tb field. Al- 
though the density fluctuations are nearly gaussian on the rel- 
evant scales, the H II regions are not. We now consider the 
pixel distribution function (PDF) in the two models; i.e. the 
fraction of spherical "pixels" of a given mass ;« p ; x that have 
brightness temperature 5Tb- We showed how to compute this 
for the inside-out model in Paper I. 7 The same technique ap- 
plies to the void model once we modify the excursion set bar- 
rier appropriately. 8 We show some example PDFs in Figure^J 
they assume m p ; x = 10 13 M or an angular resolution ~ 2.7'. 
The Figure explicitly demonstrates that although both reion- 
ization models increase the variance in 5Tb, they do so in qual- 

6 As described in Paper I, our model probably begins to break down at this 
point because the internal structure of the H II regions becomes important. 

7 Note that we do not include overlap with large bubbles in the results 
shown here, so these PDFs cannot be directly compared with observations; 
see Paper I. 

8 In this case xh (5) is a monotonic function, so the resulting PDF is non- 
singular. 




<5T b (mK) 

FIG. 4. — The PDF (top panel) and cumulative distribution function 
P(> (57),) (bottom panel) in our three models, assuming £ = 40 and m p j x = 
10 13 Mq. The solid, dashed, and dotted curves are for the inside-out model 
of Paper I, the outside-in void model, and a uniformly ionized IGM, respec- 
tively. In the bottom panel, the three sets of curves are for z = 20 (xh = 0.98), 
Z = 16 (xh = 0.88) and z = 13 (x H = 0.5), from right to left. We only show the 
PDFs for the latter two cases. 



itatively different ways. For inside-out reionization, high den- 
sity regions are highly ionized and 5Tb has a maximum value 
5T m!a on any given scale. The pixels tend to cluster strongly 
around this maximum value because 5Tb is a weak function 
of density in the rejevant range; the PDF has large variance 
because the mean 5T < 5T lmx thanks to fully ionized regions. 
For outside-in reionization, on the other hand, high-density 
pixels remain mostly neutral, while the peak and the low den- 
sity tail get stretched to small 5Tb- The result is not so far from 
a gaussian, but xh(5) has dramatically increased the variance. 
The bottom panel shows the cumulative distribution functions 
P(> 5Tb) for three redshifts, which is often useful in statisti- 
cal tests. Note that at early times, when xh is small, the dis- 
tributions closely resemble each other, but significant features 
appear even when the ionized fraction is only ~ 10%. Also, 
the inside-out model has p(5Tb = 0) ~ 0.2 at z = 13, because 
the characteristic bubble size is near m plx . Figure |4] shows the 
necessity of identifying useful tests for non-gaussianity in the 
signal. We intend to explore such tests in the future. 

5. DISCUSSION 

We have examined what 21 cm tomography of the high- 
redshift Universe can teach us about the ionization history. 
Using the reionization model described in Paper I, we pre- 
dicted the angular power spectra of the brightness temperature 
5Tb in several qualitatively different reionization scenarios. 
We have shown that 21 cm measurements offer a powerful 
probe of the reionization process, even without high signal-to- 
noise maps. Not only can they easily determine the timing of 
reionization, but particular ionization histories have measur- 
able effects on the power spectrum. For example, we showed 
that the "bubble" feature in the power spectrum disappears if a 
substantial fraction of the ionizations are due to hard photons 
that can travel large distances from their source. Also, we con- 
trasted the predictions of two different prescriptions for the 
transition between multiple generations of sources. In one, 
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the transition occurs in an approximately uniformly ionized 
medium. In the other, the transition occurs while the original 
H II regions are still growing. We showed that these differ- 
ent prescriptions lead to qualitatively different power spectra, 
even long after the first sources have shut off, because the 
first generation bubbles can imprint a persistent feature on the 
power spectrum. This could be especially important if the 
first phase of reionization ends at redshifts inaccessible to the 
next generation radio telescopes, because of either radio fre- 
quency interference or simple noise considerations (note that 
the background is approximately proportional to [1 +z] 2 7 for 
these measurements). Finally, we contrasted our "inside-out" 
reionization model with predictions for an "outside-in" model 
in which voids are ionized before dense regions. While the 
power spectrum of these models differ in detail, we showed 
that the underlying distributions of pixel brightness tempera- 
ture for a given angular resolution are dramatically different. 
We thus suggest that more sophisticated statistical measures 
of the data will be able to strongly constrain the underlying 
model of reionization. 
Of course, we have made a number of simplifying assump- 

9 See http://astrophysics.phys. emu. e du/~jbp for details on PAST. 

10 See http://www.lofar.org for details on LOFAR. 

" See http://www.skatelescope.org for details on the SKA. 



tions in all of these models. Most notably, we have ignored 
recombinations, radiative transfer, and the anisotropic distri- 
bution of sources; we discuss these issues in more detail in Pa- 
per I. Addressing these complications in a satisfactory manner 
will probably require numerical simulations. However, none 
are likely to alter our conclusions: 21 cm measurements con- 
tain an unprecedented wealth of information about the timing 
and morphology of reionization. Our simple model shows 
that the qualitative statistical features of the signal can be re- 
lated to the source and IGM properties. The next generation 
of low-frequency radio telescopes, such as the Primeval Struc- 
ture Telescope, 9 the Low Frequency Array, 10 and the Square 
Kilometer Array, 1 1 should be able to measure these statistical 
properties (ZFH04) and place strong constraints on the reion- 
ization era. 
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Science and Engineering. 
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